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Recent Developments of Domain- wall/Overlap Fermions for Lattice QCD 

Ting-Wai Chiu a * 

a Physics Department, National Taiwan University, Taipei, Taiwan 106, Taiwan. 

I review the lattice formulations of vector-like gauge theories (e.g. QCD) with domain- wall/overlap fermions, 
and discuss how to optimize the chiral symmetry for any finite N s (sites in the fifth dimension). In this formulation, 
quark propagators in gauge background can be computed efficiently through the effective 4D lattice Dirac operator. 
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1. Introduction 

The basic idea of domain-wall fermions 
(DWF) [Tl2"] is to use an infinite set of coupled 
Dirac fermion fields {ip s (x),s <E (— oo,oo)} with 
masses behaving like a step function m(s) = 
m0(s) such that Weyl fermion states can arise as 
zeromodes bound to the mass defect at s = 0. 
However, if one uses a compact set of masses, 
then the boundary conditions of the mass (step) 
function must lead to the occurrence of both left- 
handed and right-handed chiral fermion fields, 
i.e., a vector-like theory. 

Exact chiral symmetry on the lattice was pio- 
neered by Kaplan 3 with his proposal of domain- 
wall fermions on the 5-dimcnsional lattice, in 
which the fifth dimension (internal flavor space) 
is discretized with N s sites (flavors) and lattice 
spacing a§. Although the initial motivation was 
to provide a nonperturbative formulation of chi- 
ral gauge theories, the idea turns out to be natu- 
ral and well-defined for vector-like gauge theories 
(e.g. QCD), with quark fields constructed from 
the boundary modes with open boundary condi- 
tions j 1I5| . Soon after Kaplan proposed DWF for 
chiral gauge theories, Narayanan and Neuberger 
observed that the chiral determinant can be 
written as the inner-product ("overlap") of two 
fermionic many body states of two bilinear Hamil- 
tonians. For vector gauge theories like QCD, the 
fermion determinant is the product of a com- 
plex conjugate pair of chiral determinants, thus 
it is gauge invariant, real, non-negative, and the 



corresponding lattice Dirac operator for massless 
quarks (i.e., the overlap Dirac operator [7]) can be 
represented by a finite matrix of fixed shape re- 
gardless of the topology of the background gauge 
field, without undesired doubling or any fine- 
tuning. Mathematically, the overlap Dirac op- 
erator is exactly equal to the effective 4D lattice 
Dirac operator (for internal fermions dressed with 
pseudofermions) of DWF (with m q = 0) in the 
limit N s — > oo followed by a$ — > 0, 



D = m 1 + 7.5- 



H w = 75 D v 



(1) 
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where D w is the standard Wilson Dirac operator 
plus a negative parameter —too (0 < too < 2). 

For lattice QCD with DWF, in practice, one 
can only use a finite number (N s ) of lattice Dirac 
fermion fields to set up the domain wall, thus the 
chiral symmetry of the quark fields (in the mass- 
less limit) is broken. Also, the discretization in 
the fifth dimension introduces 05 into the the- 
ory. Presumably, only in the limit N s — > 00 and 
05 — > 0, the correct effective 4D theory with ex- 
act chiral symmetry can be recovered. Since N s 
and 05 are independent parameters, it might hap- 
pen that even in the limit N s — > 00, the massless 
quark propagator is chirally symmetric but still 
depends on 05, and similarly for the fermion de- 
terminant. Thus, even if 05 is an irrelevant pa- 
rameter (which decouples in the continuum limit), 
it is interesting to see whether a§ can be elim- 
inated completely at finite lattice spacing a, or 
whether it can be turned into a set of parameters 
("variable spacings") such that the chiral symme- 
try can be preserved optimally. 
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It turns out that for the conventional DWF 
with open boundary conditions 0], its quark 
propagator and effective 4D Dirac operator do not 
possess the optimal chiral symmetry for any fi- 
nite N s , and they depend on 05 even at N s — 00, 
through the Hermitian operator H in the 
transfer matrix. 

In this talk, I focus on the topics how to con- 
struct DWF action for lattice QCD |8l9ll0| such 
that the quark propagator and the effective 4D 
lattice Dirac operator possess optimal chiral sym- 
metry for any finite N s and gauge background. 
Due to the time constraint, I could not review 
recent numerical results on lattice QCD with 
domain- wall/overlap fermions. 

2. Lattice QCD with conventional DWF 

First, we examine the action 2 of conventional 
DWF with open boundary conditions | II oj 

N s 

s,s'=l x,x' 

+5 x , x ,D 5 (s,s')}iP(x',s'), (2) 

where D w is the Wilson-Dirac operator minus a 
parameter mo (0 < 057710 < 2, < mo < 2), and 



D 5 (s,s') = S s 



P_<5 



P+5 



+ O s ',s-l, 



P± = -(l±7s), 

with boundary conditions: 

P+^(x,0) = -rm q P+ip(x,N s ), (3) 
P_ip(x, N s + 1) = -rm q P-i){x, 1), (4) 
r = [m (2 - B 5 m )p 1 . (5) 

Here m q is the bare quark mass, and the quark 
helds coupling to physical hadrons are con- 
structed from the boundary modes: 

q(x) = Vf[P-ip(x, 1) + P+fl>(x, N.)] (6) 
q(x) = SF${x, l)P+ + i>(x, N S )P_}. (7) 

To regularize the fcrmion determinant, Pauli- 
Villars (pseudofermion) fields {4>xs,4'xs} carrying 

2 Here we suppress the lattice spacing a, and the Dirac, 
flavor, and color indices. 



the same (flavor, color, Dirac) indices of {ip,4>} 
but obeying the Bose statistics are introduced [111 
[H] . The action of pseudofermions Ap f is exactly 
the same as @ except replacing {tp, by {</>, </>}, 
with boundary conditions (rm q = 1): 

P+0(x,O) = -P+<f>(x,N s ), 
P-4>(x,N s + l) = -P_(j,(xA). 

Thus the total action of the system is 

A = A g + Af + A pf 

where Ag is the gluon action. 

Then the fermion determinant, the effective 4D 
lattice Dirac operator J21> and the quark prop- 
agator ^3] m background gauge field can be ob- 
tained by integrating over all heavy fermion fields 
and pseudofermion fields. In general, the gener- 
ating functional for n-point Green's function of 
quark fields q and q can be derived as |10| 

J[dU]e- AG W det£»(m g ) e J ( fl ^')" lj 



Z[J, J] = 



J[dU]e- A °M detD(m q 



where J(x) and J(x) are the Grassman sources 
of q(x) and q(x) respectively, and 

1 + lsS(H) 



rD c = 
S(H) 

T - 

H 

D(m q )- 



. J<N. 



1 + T N > 
1 - a 5 H 



a 5 HR(a 2 5 H 2 



1 

75 



a 5 H 
D m 



a$D w 
(D c + m q ){l 4 
: [1 + rm q + (1 



(8) 
(9) 
(10) 

(11) 

rDA- 1 

-rm,)7isS(£0]/2r(12) 



D{m q ) 1 = (l-rm q )(D c + m q ) l +r 



(13) 



The most remarkable feature of DWF is the 
emergence of the quark propagator (D c + m 9 ) _1 
which preserves all vital symmetries of its coun- 
terpart [7^(9^1 + igA^) + niq]^ 1 in continuum, 
while the sea quark propagator D{m q )~ 1 for in- 
ternal quark loops is equal to (D c + rriq)^ 1 times 
a constant factor (1 — ram q ) plus a constant ra, 
as shown in Eq. Ijl^JI . 

In the limit N s -> 00, S(H) -> ^L= = sgn(P), 
then the quark propagator (D c + m g ) _1 in the 
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massless limit (m q — > 0) is chirally symmetric. 
Consequently, the effective 4D lattice Dirac op- 
erator for massless internal quark loops, -D(O) = 
D c (l + rDc)- 1 = [1 + lb S{H)]/2r = D, satisfies 
the Ginsparg- Wilson relation ^1] 

Dj 5 + l5 D = 2rD l5 D. 

Further, taking 05 — > 0, then S(H) — ► sga(H w ), 
and D is exactly equal to the overlap Dirac oper- 
ator ft), (But a 5 mwsi fee nonzero at finite N s ). 
Several remarks are as follows. 

(i) The axial anomaly is recovered through the 
Chern-Simmons current |15I4| . 

(ii) In the limit N s — > 00, the non-singlet flavor 
symmetry SU^Nf) x SUn(Nf) of 2V/ massless 
quarks is exact at finite lattice spacing 0j. 

(iii) The effective 4D lattice Dirac opera- 
tor D(m q ) is exponentially-local for sufficiently 
smooth background gauge fields |16I17I18) . 

(iv) Any quark observable in QCD can be ob- 
tained from Z[J, J] by differentiation, and it pos- 
sesses the discrete symmetries (C, T, P) of its 
counterpart in continuum. 

(a) Quark propagator 



S 2 Z[J,J] 



-A G 



n 



5J(x)5J(y) 
det£>(m/) (D 



j=j=o 
. + m 



1 

<llx.y 



f[dU}e-^H fdet D(m f ) 
(b) Current-current correlator 



(d(x)~f 4 P-s(x)s(0)j 4 P-d(0)} 



:1aP- 



5J d {x)'* ~ 5J s {x) SJ s {0) ~ SJ d (0) 
J[dU}e- Ao Uf deW(m f )0 K (x) 
J[dU]e~ A ol\ f detD(m f ) 



Z[J, J] 



where 

K {x) = -tr[(D c +m d ) -^7 4 P_( J D c -|-m 5 .)- 1 74P_]. 

Note that the V-A structure of the left-handed 
quark currents is preserved exactly. 

(v) If one uses the effective 4D lattice Dirac 
operator D(m q ) (or any GW Dirac operator) to 
define the quark action, then subtleties might 
emerge. (Recall that CP in chiral gauge theories 
with GW fermion is broken by 0{a) effect 



Consider the massless (m q — 0) case, 
A F = ^2 q(x)D x , y q(y) = qDq 

=J2[q(x)P+D x , y P-q(y) + q{x)P-D x , y P + q{y)] 

where P± = |[1 ± 75(1 - 2rD)} and P± = (1 ± 
75)/2 are the chiral projectors for the quark fields 
q and q respectively. Then the left-handed quark 
current does not manifest the V-A structure, e.g. 
d(x)j fi P-s(x) ^ d(x)P+ 7^5(0;). Consequently, 
the effective weak Hamiltonian constructed from 
these left-handed quark currents would break the 
SUl(2) gauge symmetry by 0(a) effect [2D] . 

Note that this problem cannot be resolved by 
re-defining the quark field as q = (1 — rD)q = 
(1 + rD c )~ 1 q, since Af = qDq — qD c q, where D c 
becomes nonlocal as N s — > 00. 

Thus, for any Ginsparg- Wilson fermion action, 
it is impossible to define local quark fields with 
exponentially-local kernel such that the correla- 
tion function of quark fields in gauge background 
can be expressed as 

(q(x)q(y)) = (D c + m,)^ . 

Even though (D c + m q )~ l can be computed via 
D{m q )~ 1 , conceptually, the valence quark fields 
have to be defined in terms of the boundary modes 
of domain-wall fermions. 

Nevertheless, the conventional DWF has its 
own deficiencies: 

(a) For any finite N s , S(H) is not the optimal 
rational approximation for sgn(77) 8], thus the 
chiral symmetry is not preserved optimally. 

(b) Even for N s = 00, S(H) = H/VW still 
depends on 05 through H = ^D w {2 + a^Dw)^ 1 . 

In the following, I discuss how to solve these 
two problems. 

3. Zolotarev optimal rational polynomial 

The deviation of S(H) from sgn(JJ) can be 
measured in terms of 



o~(S) = max 



Y^{sgn(H) - S{H)}Y 



< max |sgn(ry) 



yty 
s(v)\ 



(14) 
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where {77} are eigenvalues of a§H . Using the sim- 
ple identity 

|sgn(ar) - S(x)\ = |1 - V^R(x 2 )\ 

where S(x) = xR(x 2 ), we can rewrite (|14|) as 



and 



a(S) < max 
" {v 2 } 



1 



(15) 



where {r/ 2 } are eigenvalues of a\H 2 . Clearly 
the problem of finding the optimal rational 
approximation S op t(x) of sgn(x) with x G 
[ %maxi Xmin] U [xmin, x raax \ is equivalent 
to finding the optimal rational approximation 
R op t{x 2 ) ofix 2 )- 1 / 2 withx 2 £ [x 2 min ,x 2 max ]. 

According to de la Vallee-Poussin's theorem 
and Chebycheff's theorem. 2], the necessary and 
sufficient condition for an irreducible rational 
polynomial (where m > n, Pi,qi > 0) 

(n vri) [ \ PnX n + Pii-lX^ 1 -\ h PO 



q m x m + q m -ix m 1 H h q Q 



to be the optimal rational polynomial of the in- 
verse square root function x -1 / 2 , < x m i n < 
x < x m ax is that S(x) = 1 — \fx A n ' m \x) has 
n + r7i + 2 alternate change of sign in the interval 
[x m im x m ax], and attains its maxima and minima 
(all with equal magnitude), say, 

S(x) = -A, +A, ■ • • , (_i)«+™+2a 

at consecutive points (a;,, i = 1, • • • , n + m + 2) 

Xmin X\ <C X2 \ ' • • <C — Xmax • 

In other words, if r^™'™) satisfies the above con- 
dition, then its error 



(,■) 



.( r (i»,m)j = max l-^fM 

XG [imiTi^mai] 

is the minimum among all irreducible rational 
polynomials of degree (n, m). 

It has been shown [8] that i?(x 2 ) JJjJ of the 
conventional DWF (i.e., the polar approximation) 
does not satisfy above criterion, thus is not the 
optimal rational approximation for (x 2 ) -1 / 2 . 

The optimal rational approximation for inverse 
square root function (x 2 ) _1 / 2 ,x 2 S [1,6] was first 
obtained by Zolotarev in 1877, using Jacobian el- 
liptic functions. Explicitly, 



n , 

R ( r\x 2 )=d n \\-±_ 

1=1 



X 2 /C 2 l 



+ X 2 /c 2 /_l 



(16) 



R 



(n— l,n) / 2 



(x 2 )=<z" 



n;L 1 1 (i+.x 2 /^) 
nr=i(i+^ 2 /^-i) 



(17) 



where the coefficients do, d , ci and c\ are ex- 
pressed in terms of elliptic functions |21| with ar- 
guments depending only on n and b. A detailed 
discussion of Zolotarev's result can be found in 
Akhiezer's two books [22- ( see also a recent dis- 
cussion |22] in the context of lattice QCD.) 

The first application of Zolotarev optimal ra- 
tional polynomial in lattice QCD was to approx- 
imate [231 (-^ui) -1 ^ 2 m f ne overlap Dirac oper- 
ator which improves the (optimal) rational 
approximation |24j via the Remes algorithm 3 . 

For lattice QCD with DWF at finite N s , the 
optimization problem is how to construct a DWF 
action such that the operator S (H) © is replaced 
with 



S op t{H) 



r hR%' n \h 2 ), 



N s = 277 + 1, 

N s = 277, 



(18) 



where h = H/X min , and b = A^ oa ./A^ in (A min 
and X max are the minimum and the maximum of 
the eigenvalues of \H\). However, this could not 
be solved for the conventional DWF, due to the 
functional form of H = jsD w (2 + a^Dw)^ 1 . 

4. Borici's variant of DWF 

For the effective 4D Dirac operator D(m q ) (|12|) 
in the fermion determinant, it turned out to be 
possible to modify the conventional DWF action 
such that its S{H) is replaced by S(H W ). 

The prescription |25j is to replace S XtX iD 5 (s, s') 
in J2J with 



{a 5 D w - l) x ,x'(P-6s',8+i + P+S S ',s-i). 



The boundary conditions ©-Q) and quark fields 
©-0 are the same except replacing r © with 
r = 1/(2too). After introducing pseudofermion 

3 The Remes algorithm is a numerical (iterative) scheme to 
implement the criterion of de la Vallee-Poussin's theorem 
and Chebycheff's theorem to obtain optimal rational ap- 
proximation of any continuous function f(x). However, in 
the case (a; 2 ) -1 ' 2 , Remes algorithm must give less precise 
numerical coefficients than those computed directly from 
Zolotarev's exact solution. 
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fields (rm q — 1), the fermion determinant can be 
evaluated as 

det{T> F V PF ) = detD 
= det{m g + (m - m q /2)[l + l5 S(H w )}}. (19) 

In the limit N s oo, S(H W ) -> H w {Hl)~ x / 2 , 
then the fermion determinant l|19l) is independent 
of 05. 

However, unlike the conventional DWF, its 
quark propagator in gauge background is ab- 
normal, 

(q(x)q(y))^(D c + m q )-i, (20) 

thus its quark propagator as well as other quark 
observables do not manifest the discrete symme- 
tries of their counterparts in continuum, simi- 
lar to the problem encountered in the Ginsparg- 
Wilson fermion. 

5. Optimal domain-wall fermion (ODWF) 

The problem (|20[1 can be solved JUj by ap- 
pending two new boundary layers at s — and 
s = N s + 1, both with the constraint a§ — 0, and 
re-defining the quark fields with these new bound- 
ary modes. 

Further, one can turn 05 at each layer into a 
parameter u s . Then this set of parameters {uj s } 
can be constructed such that S(H W ) is equal to 
S pt{H w ), the optimal rational approximation for 
sga(H w ) 0. 

The action of optimal domain- wall fermion |1U| 
can be written as 

N s + 1 

A F = ^V odw ^ ^^ XS {{^ S D W + 1) XX >5 SS , 

s,s'—0 x,x' 

+ (uj s D w - l) xx ,(P + 8 s ^ s -i + P_6 S ', s +i)}ipx's'(21) 
with boundary conditions 

P+i/)(x,-i) = -rm q P + ip(x,N 3 + l), 
P_ip(x,N s + 2) = -rm q P-<tp(x,0), r = 



1 



2mo 

where ujq — wjv s +i — 0, and other {co s } are fixed 
as follows. The quark fields are constructed as 

q(x) = Vr~[P-^(x,0)+P+^(x,N s + l)} 
q{x) = Vr~$(x,0)P++^(x,N s + l)P^} . 



Then the quark propagator in gauge background 
can be derived as ITUl 



(q(x)q(y)} = (D c + m,)"* 
where 

l+ l5 S (H w ) 



D c = 2m 



1 - j!yS (H w ) ' 



l-Y\ Ns T 

Q / If \ 1 Ls=l » rp 

1 + lls=l ls 



(22) 

(23) 
(24) 



Now requiring S (H W ) equal to the optimal ratio- 
nal approximation of sgn(H w ) amounts to solving 
the parameters {lo s } from the equation 



{, n(n,n)/,2\ 
h w R^ n \hl), N s = 2n. 



2n+ 1, 



This is equivalent to solving the roots (u s = 
ujJ 2 , s = 1, • • • ,N S ) from the nonlinear equations 



0, 



1- v^4 n ' n) o 

l-VSH?~ 1,n) («)=0, N s 



2n 
2n. 



(25) 



The exact solution of 112311 has been obtained in 
|U, which gives 



■y/i -«' 2 sn 2 (v s ;k'), s = l,.-.,JV,(26) 



Here sxi(v s ;k') is the Jacobian elliptic function 
with modulus k' = \fl — A 2 „ in /A 2 la:I ,, and argu- 
ment u„ 



where 



m = n 



1 


^3A 




hA)3 







A' + 



2K' 



n ° 2 (^ 



1=1 sn- 

N. 



2 ( (2l-l)K> , 



A' = Vl-A 2 , (27) 



if' is the complete elliptic integral of the first kind 
with modulus k', and is the elliptic theta func- 
tion. From (PI), one has X m \ x < uj s < A~* n , 
since sn 2 (;) < 1. Note that A m i n (A max ) can be 
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fixed to be the greatest lower bound (least upper 
bound) of the eigenvalues of \H W \ for the set of 
gauge configurations under investigation. Clearly 
the ODWF action ^21)) is always ultralocal, no 
matter how one fixes X m in o,nd \ ma x- 

It should be pointed out that the new bound- 
ary layers at s = and s = N s + 1 with ujq = 
ujn 3 +i = were first introduced in JOjj several 
months after the weights H26|) had been obtained 
in [H| . The purpose of these two new boundary 
layers is to render the quark propagator satisfying 
()22l) such that any quark observable manifests the 
discrete symmetries of its counterpart in contin- 
uum. Since To = Tjv s +i = 1, the transfer matrix 
as well as S (H W ) are not affected by these two 
new boundary layers. 

After introducing the pseudofermion fields 
(m q = 2mo), the generating function for n-point 
function of q and q can be derived as jlOj , 

- _ f[dU}e- A °M detD(m 9 ) e J (^+ m ^ lj 
Z[J ' JJ ~ J[dU}e--*G[u] detD(m q ) 

where 

JJ C = zmo - 

D(m q )=m q + (m - m q /2)[l + j 5 S op t{H w )} (28) 

It has been shown that D(m q ) (for any m q and 
N s ) is exponentially- local for sufficiently smooth 
gauge backgrounds UJ. Note that for any finite 
N s , D(0) (|28|l is exactly equal to the overlap 
Dirac operator with (if^,) -1 / 2 approximated by 
Zolotarev optimal rational polynomial, (|16|l (N s 
= odd), or jni> {N a = even). 

Since S pt(H w ) is the optimal rational approx- 
imation for sgn(iJ lu ), its error is bounded uni- 
formly over the entire range of [1, A^^/A^J, 
and is the minimum among all rational polyno- 
mials of the same degree, 

<j(S opt ) <\^^ A(b)e- C ^ N % b = X 2 max /X 2 min 

where A is defined in 127|) . and 

A(b) ~ 4.06(l)6- aoo91 ( 1 Hn(6) a0042 ^ 
c{b) ~ 4.27(45)ln(5)- a746 ( 5 ) 

which can be estimated by asymptotic expansion 
and numerical evaluation of A. 



Thus for any set of gauge configurations, one 
can determine what values of N s and b (i.e., how 
many low-lying eigenmodes of should be pro- 
jected out) are required to attain one's desired 
accuracy in preserving the chiral symmetry. 

The quark propagator in gauge background 

{q(x)q(v)) = (d c + m q )-i 

=J2(P-5 s .o + P + S s>Na+1 )V-^ w (x, s; y, s') x 

s,s f 

(P+tf^O + P-^AT.+l) 

can be computed in two different ways: 

(a) To solve the linear system of the 5D lattice 
Dirac operator 

T>odw{x, s; y, s')Y = (P+S sfi + P-5 SiNs+ i) ■ I 

(b) To solve the system of the effective 4D lattice 
Dirac operator 

D{m q )D\m q )Y = 1 
= {m 2 + (m 2 - m 2 /4) [2 + ( 75 ± l)S opt (H w )]}Y 

with conjugate gradient, where the matrix-vector 
product S op t(H w ) - Y can be obtained by invoking 
another conjugate gradient (inner CG). It turns 
out that the best algorithm to compute S op t ■ Y 
is Neuberger's double pass algorithm |27j. 

At present, the scheme (b) seems to be much 
more efficient than the scheme (a), in terms of 
the CPU time and the memory space. Further, 
there are multi-shift CG algorithms available for 
the scheme (b) to compute quark propagators for 
a set of quark masses, but so far there is no multi- 
mass algorithms for the scheme (a). For a recent 
quenched QCD calculation via the scheme (b), 
see |25) . 

At this point, it is instructive to point out the 
remarkable features of Neuberger's double pass: 

(i) the memory storage for the conjugate gra- 
dient is constant (5 vectors), independent of the 
degree n of the rational polynomial i?( n_1 >™). 

(ii) the CPU time is almost independent of n, 
thus the matrix-vector product R ■ Y can be ap- 
proximated to very high precision (at large n) 
without noticably extra costs. 

(iii) there exists a threshold nx such that the 
double pass algorithm is faster than the single 
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pass algorithm for n > rij- (where nj* — 12 — 25 
for most platforms). 

Note that the last two properties have just been 
unveiled recently (2H|. Undoubtedly, Neuberger's 
double pass algorithm will become the standard 
for computing R(H 2 ) ■ Y in lattice QCD. 

6. Conclusions 

A new DWF action (J2U for lattice QCD has 
been constructed such that the quark propaga- 
tor and the effective 4D Dirac operator for inter- 
nal quark loops have optimal chiral symmetry for 
any N s and gauge background. The quark fields 
are constructed from the new boundary modes at 
s = and s — N s + 1 with u>q = lon 3 +i = such 
that the quark propagator satisfies (|22|l and any 
quark observable manifests the discrete symme- 
tries of its counterpart in continuum. The quark 
propagator in gauge background can be computed 
efficiently through the effective 4D lattice Dirac 
operator, and the dynamical quarks can be sim- 
ulated similar to the Wilson fermion in 5D. 
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